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Abstract: In this article we integrate in closed implicit form the Navier- 
Stoker equations for an incompressible fluid in a smooth compact manifold with- 
out boundary, and in particular, in a compact manifold which is isometrically 
embedded in Euclidean space, and finally in Euclidean space itself. We further 
integrate the kinematic dynamo problem of magnetohydrodynamics, i.e. the 
equations of passive transport of a magnetic field on a fluid. We carry out these 
integrations through the application of the methods of Stochastic Differential Ge- 
ometry, i.e. the gauge theoretical formulation of diffusion processes on smooth 
manifolds. Thus we start by defining the invariant infinitesimal generators of dif- 
fusion processes of differential forms on smooth compact manifolds, in terms of 
the laplacians (on differential forms) associated with the Riemann-Cartan-Weyl 
(RCW) metric compatible connections. These geometries have a torsion ten- 
sor which reduces to a trace 1-form, whose conjugate vector field is the drift of 
the diffusion of scalar fields. We construct the diffusion processes of differential 
forms associated with these laplacians by using the property that the solution 
flow of the stochastic differential equation corresponding to the scalar diffusion 
is -assuming Hoelder or Sobolev regularity conditions- a random diffeomorphism 
of the manifold. We apply these constructions to give a new characterization of 
the Navier-Stokes equation for the velocity one-form of an incompressible fluid 
as a non-linear diffusion process determined by a RCW connection. We prove 
the Navier-Stokes equations to be equivalent to a linear diffusion equation for the 
vorticity and the Poisson-de Rham equation for the velocity with the vorticity 
as a source, extending thus to manifolds, the approach in practice in the random 
vortex method in Computational Fluid Dynamics. We give the invariant random 
stochastic differential equations for the position (as a Lagrangian representation) 
of the fluid particles and thus obtain a random diffeomorphism which is a solu- 
tion of the Navier-Stokes equation. We solve the Cauchy problem for the heat 
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equation for the vorticity two-form and the Dirichlet problem for the Poisson-de 
Rham equation for the velocity one-form in two different instances. Firstly, using 
the gradient flows obtained by a Nash isometric embedding of the manifold in Eu- 
clidean space, and secondly, for the more general case of an arbitrary manifold, by 
running curvature and fluid-deformation tensor dependant random flows on the 
tangent bundle, projecting on the random fluid' particles flow. We implement our 
general construction on Euclidean space by simply taking the former construc- 
tion under the identity embedding. We extend these methods to the integration 
of the kinematic dynamo problem. From the fact that any metric compatible 
connection can be constructed in terms of the push-forward connection and us- 
ing simple stochastic analysis considerations, we prove that any diffusion process 
generated by a RCW connection admits a random flow representable through 
a purely diffusive process (i.e. zero drift), in any dimension other than 1. We 
apply this to prove that the random flow of a viscous fluid obeying the classical 
Navier-Stokes equations on a smooth manifold of dimension other than 1, can be 
represented as a purely diffusive process, where the new diffusion tensor can be 
constructed in terms of the velocity and the original diffusion tensor associated 
to the metric and the kinematical viscosity. A similar construction we prove to 
be valid for the fluid flow of a passively magnetized viscous or inviscid fluid. 
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1 Introduction. 



The purpose of this article is two-fold: Firstly to give implicit random rep- 
resentations for the solutions of the Navier-Stokes equation for an incompressible 
fluid and for the kinematic dynamo problem of magnetohydrodynamics, in several 
instances; firstly, on an arbitrary compact orientable smooth manifold (without 
boundary) -following our presentation in [58,59]- and further in the case in which 
the manifold is isometrically embedded in Euclidean space, to finally give the 
expressions for Euclidean space itself. Secondly, to present as a basis for such 
an integration, the gauge-theoretical structures of Brownian motion theory and 
the stochastic analysis rules associated to them. Thus, the method of integra- 
tion we shall apply for our objectivs stems from stochastic differential geometry, 
i.e. the gauge theory of Brownian processes in smooth manifolds and Euclidean 
space developed in the pioneering works by Ito [15], Elworthy and Eells [13], P. 
Malliavin [11] and further elaborated by Elworthy [12], Ikeda and Watanabe [14], 
P. Meyer [40], and Rogers and Williams [37]. Associated to these geometrical 
structures which can be written in terms of the Cartan calculus on manifolds 
of classical differential geometry, we shall present the rules of stochastic analysis 
which describe the transformation of differential forms along the paths of gener- 
alized Brownian motion generated by these geometries, setting thus the method 
for the integration of linear evolution equations for differential forms; this is the 
well known martingale problem approach to the solution of partial differential 
equations on manifolds [30]. 

While classical Hamiltonean systems with finite degrees of freedom may ap- 
pear to have a random behavior, in fluid dynamics it is known that the Euler 
equation for an inviscid fluid is a hamiltonean system with infinite degrees of 
freedom supporting as well infinite conserved quantities; such a system appears 
to be non-random [56] . The situation is radically different in the case of a viscous 
fluid described by the Navier-Stokes equations. In this case, there is a second- 
order partial derivative associated to the kinematical viscosity, which points out 
to the fact that there is a diffusion term, which can be thougth as related to a 
Brownian motion. Thus in the viscous case, there is from the very beginning 
a random element. While in the Euler case the group of interest is the (Rie- 
mannian) volume preserving of diffeomorphims, it will turn out in the course of 
these studies, that there is an active group of random diffeomorphims which 
represent the Lagrangian random flow of the viscous fluidparticles. In this case, 
when there is a non-constant diffusion tensor describing the local amplification of 
noise, these diffeomorphisms do not preserve the volume measure, contributing 
at a dynamical level -as it will turn out- to the complicated topology of turbulent 
and magnetized flows [56]. 
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The essential role of randomness in Fluid Dynamics appears already at an 
experimental level. The analysis of the velocimetry signal of a turbulent fluid 
shows that its velocity is a random variable, even though that the dynamics is 
ruled by the deterministic Navier-Stokes equation [8]. The concept of a turbulent 
fluid as a stochastic process was first proposed by Reynolds [16], who decomposed 
the velocity into mean velocity plus fluctuations. The Reynolds approach is cur- 
rently used in most numerical simulations of turbulent fluids in spite of the fact 
that it leads to unsurmountable non-closure problems of the transport equations.; 
see Lumley [18], Mollo and Christiansen [38]. Furthermore, the Reynolds decom- 
position is non invariant alike the usual decomposition into drift and white noise 
perturbation in the non-invariant theory of diffusion processes. Other treatments 
of stochasticity in turbulence were advanced from the point of view of Feynman 
path integrals, as initiated by Monin and Yaglom [17]. From the point of view of 
diffusion processes, invariant measures for stochastic modifications of the Navier- 
Stokes equations on euclidean domains, have been constructed by Vishik and 
Fursikov [36] and Cruzeiro and Albeverio [42]. (It is important to remark that 
the existance of an invariant measure for NS as a classical dynamical system is 
the starting point of the classical dynamical systems approach to turbulence; see 
Ruelle [48].) Contemporary investigations develop the relations between random- 
ness and the many-scale structure of turbulence which stems from the Kolmogorov 
theory as presented by Fritsch[8] and Lesieur [3], and apply the renormalization 
group method; see Orzsag [19]. 

A completely new line of research followed from the understanding of the 
fundamental importance of the vorticity (already stressed by Leonardo da Vinci) 
in the self-organization of turbulent fluids, which was assessed by numerical sim- 
ulations by Lesieur [3,4], and theoretically by Majda [39] and Chorin [1]. It 
was found that the Navier-Stokes equations for an incompressible fluids on Eu- 
clidean domains yields a linear diffusion equation for the vorticity which becomes 
a source for the velocity through the Poisson equation: Solving the latter equa- 
tion, we can obtain an expression for the fluids velocity in 2D. This observation 
was the starting point for the random vortex method in Computational Fluid 
Mechanics largely due to Chorin [1,2,6]. This conception lead to apply methods 
of statistical mechanics (as originally proposed by Onsager [20]) to study the 
complex topology of vortex dynamics and to relate this to polymer dynamics 
[1]. In the random vortex method a random lagrangian representation for the 
position of the incompressible fluid particles was proposed. Consequently, the 
Navier-Stokes linear (heat) equation for the vorticity was integrated only for two 
dimensional fluids (implicit to this is the martingale problem approach quoted 
above), while the general case was numerically integrated by discretization of the 
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this heat equation; see Chorin [1,2]. The difficulty for the exact integration in 
the general case apparently stems from the fact that while in dimension 2 the 
vorticity 2-form can be identified with a density and then the integration of the 
Navier-Stokes equation for the vorticity follows from the application of the well 
known Ito formula for scalar fields, in the case of higher dimension this iden- 
tification is no longer valid and an Ito formula for 2-forms is required to carry 
out the integration. This formula became only recently available in the works by 
Elworthy [27] and Kunita [24], in the context of the theory of random flows on 
smooth manifolds. 

The importance of a Stochastic Differential Geometry treatment of the Navier- 
Stokes equation on a smooth n-manifold M stems from several fundamental facts 
which are keenly interwoven. For a start, it provides an intrinsic geometrical 
characterization of diffusion processes of differential forms which follows from 
the characterizations of the laplacians associated to non-Riemannian geometries 
with torsion of the trace type, as the infinitesimal generators of the diffusions. 
In particular, this will allow to obtain a new way of writing the Navier-Stokes 
equation for an incompressible fluid in terms of these laplacians acting on dif- 
ferential one- forms (velocities) and two- forms (vorticities). Furthermore, these 
diffusion processes of differential forms, are constructed starting from the scalar 
diffusion process which under Hoelder continuity or Sobolev regularity condi- 
tions, yields a time-dependant random diffeomorphism of M which will represent 
the Lagrangian trajectories for the fluid particles position. This diffeomorphic 
property will allow us to use the Ito formula for differential forms (following 
the presentation due to Elworthy) as the key instrument for the integration of 
the Navier-Stokes equation (NS for short, in the following) for an incompressible 
fluid. Thus, it is the knowledge of the rules of stochastic analysis what sets the 
martingale problem approach to the solution of NS when he have transformed it 
to an equivalent system which is essentially linear. This transformation of non- 
linear equations to linear ones appears as well in an identical gauge-theoretical 
approach to quantum mechanics as a theory of Brownian motion. In this case, the 
massive Dirac-Hestenes equation for a Dirac-Hestenes operator field (conceived 
as a geometrical matter field) can be transformed in a Clifford bundle approach 
to the sourceless Maxwell equation [26]. We must remark that the geometrical 
structures on which the gauge-theoretical foundations of Brownian motion are 
introduced, belong in fact originally to the gauge theory of gravitation, including 
not only translational degrees of freedom, but additionally spinor fields [46,57]. 

For the benefit of completeness, we have included in this article an Appendix 
which reviews briefly the main probabilistic and analytical notions appearing in 
this article. This article represents the contents of a lecture delivered in December 
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18, 2000, at the Seminar in Geometry and Topology of Fluid Flows, Isaac Newton 
Institute for Mathematical Sciences, and expands on previous results [59,62,63]. 

2 Riemann-Cartan-Weyl Geometry of Diffusions 

The objective of this section is to show that the invariant definition of a 
"heat" (Fokker-Planck) operator requires the introduction of linear connections 
of a certain type. 

Let us consider for a start, a smooth n-dimensional compact orientable man- 
ifold M (without boundary), on which we shall consider a second-order smooth 
differential operator L. On a local coordinate system, (x a ),a = 1, ...,n, L is 
written as 

L = \g^{x)d a d p + B a (x)d a + c(x). (1) 

From now on, we shall fix this coordinate system, and all local expressions shall 
be written in it. 

Although formally, there is no restriction as to the nature of M, we are 
really thinking on a n-dimensional space (or space-time) manifold, and not in a 
phase-space manifold of a dynamical system. 

The principal symbol a of L, is the section of the bundle of real bilinear 
symmetric maps on T*M, defined as follows: for x € M, pi € T*M, take C 2 
functions, fi : M — > R with fi(x) = and dfi(x) = pi, i = 1, 2; then, 

a(x)( Pl , P2 ) = L(f 1 f 2 )(x). (2) 

Note that a is well defined, i.e., it is independent of the choice of the functions 
/<,» = 1,2. 

If L is locally as in (1), then a is locally represented by the matrix (g a ^). 
We can also view a as a section of the bundle of linear maps L(T*M, TM), or as 
a section of the bundle TM TM, or still as a bundle morphism from T*M to 
TM. If a is a bundle isomorphism, it induces a Riemannian structure g on M, 
g: M -> L(TM, TM): 

g(x)(v 1 ,v 2 ) :=< o-(x)~ 1 v 1 ,V2 > x , 

for x <G M, v\,V2 G T*M. Here, < .,. > x denotes the natural duality between 
T*M and T X M. Locally, g(x) is represented by the matrix \{g a ^ (x))" 1 . Consider 
the quadratic forms over M associated to L, defined as 

Qx(Px) = 7i< Px,0- x (p x ) > x , 
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for x € M, p x e T*M. Then, with the local representation (1) for L, Q x is 
represented as \{g a ^ (x)). Then, L is an elliptic (semi-elliptic) operator whenever 
for all x € M, Q x is positive-definite (non-negative definite). We shall assume in 
the following that L is an elliptic operator. In this case, a is a bundle isomorphism 
and the metric g is actually a Riemannian metric. Notice, as well, that a{df) = 
grad /, for any / : M — > R of class C 2 , where grad denotes the Riemannian 
gradient. 

We wish to give an intrinsic description of L, i.e. a description independent 
of the local coordinate system. This is the essential prerequisite of covariance. 

For this, we shall introduce for the general setting of the discussion, an 
arbitrary connection on M, whose covariant derivative we shall denote as V. We 
remark here that V will not be the Levi-Civita connection associated to g; we 
shall precise this below. Let <r(V) denote the second-order part of L, and let us 
denote by Xo(V) the vector field on M given by the first-order part of L. Finally, 
the zero-th order part of L is given by L(l), where 1 denotes the constant function 
on M equal to 1. We shall assume in the following, that L(l) vanishes identically. 

Then, for / : M — > R of class C 2 , we have 

a(V)(x) = l -trace(V 2 f){x) = \(Vdf){x)\ (3) 

where the trace is taken in terms of g, and Vdf is thought as a section of 
L(T*M, T*M). Also, X (V) = L - cr(V). If is the local representation 
for the Christoffel symbols of the connection, i.e. V = r^gfir, then the 
local representation of <r(V) is: 

a(V)(x) = l -g^{x){d a d p + iZp(x)dJ, (4) 

and 

X (V)(x) = B a (x)d a - l -g^{x)fl^. (5) 

If in particular, we take V the Levi-Civita connection associated to g, which we 
shall denote as V s , then for any / : M — > R of class C 2 : 

o(V 9 )(df) = itrace((V^) 2 /) = ^trace(W/) = -l/2div g grad / = l/2A fl /. (6) 

Here, A g is the Levi-Civita laplacian operator on functions; locally it is written 

as 

A g = g- 1 / 2 d a {(g 1 ^d p ) ] g = det(g a p), (7) 
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and divg is the Riemannian divergence operator on vector fields X = X a (x)d a : 

div g (X) = -g- l l 2 d a {g l ' 2 X a ). (8) 

Note the relation we already have used in eqt. (6) and will be used repeatedly; 
namely: 

divg(X) = -5X, (9) 

where 5 is the co-differential operator (see (23) below), and X is the one- form 
conjugate to the vector field X, i.e. X a = g a pX^. 

We now take V to be a Cartan connection with torsion [10,46], which we ad- 
ditionally assume to be compatible with g, i.e. Vg = 0. Then <r(V) = ^trace(V 2 ). 
Let us compute this. Denote the Christoffel coefficients of V as Fg 7 ; then, 

^7 = 1^} +l/2*&, (10) 

where the first term in (10) stands for the Christoffel Levi-Civita coefficients of 
the metric g, and 

K% = Tp-y + 8%, + S^p, (11) 

is the cotorsion tensor, with Sp = g au gp K T^, and Tp = Tp*^ — T^p the skew- 
symmetric torsion tensor. 

Let us consider the Laplacian operator associated with this Cartan connec- 
tion, defined -in extending the usual definition- by 

H(V) = l/2traceV 2 = l/2g af3 V a Vp; (12) 

then, cr(V) = H(V). A straightforward computation shows that that H(V) only 
depends in the trace of the torsion tensor and g: 

H(V) = l/2A g + g al3 Qpd a = H (g, Q), (13) 

with Q = Tj^pdxP, the trace-torsion one- form. 

Therefore, for the Cartan connection V defined in (10), we have that 

a(V) = itrace(V 2 ) = iA, + Q, (14) 

with Q the vector-field conjugate to the 1-form Q: Q(f) =< Q,grad / >, / : 
M — > R. In local coordinates, 

Q a = g aP Qp- 
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We further have: 



X (V) = B-I^|J U-Q, (15) 



Therefore, the invariant decomposition of L is 



^trace(V 2 ) + X (V) = ^A g + b, (16) 



with 



^"H^K (17) 

Notice that (15) can be thought as arising from a gauge transformation: b — > 6— Q, 
with 6 the 1-form conjugate to b. 

If we take for a start V with Christoffel symbols of the form 

r% = {^} + ^ Q, - 9m Q a } (is) 

with 

Q = b, i.e. Q = b, (19) 
we have in writing now V for the covariant derivative of (18) 

X (V) = 0, 

and 

H (g, Q) = <r(V) = ^trace(V 2 ) = itrace((V^) 2 ) + Q = X - A g + b. (20) 

Therefore, for V as in (18) we obtain a gauge theoretical invariant representation 
for L given by 

L = H(V) = l -V 2 = (t(V) = ^trace((V^) 2 ) + Q = H (g, Q). (21) 

The restriction we have placed in the metric-compatible V to be as in (18), 
i.e. only the trace component of the irreducible decomposition of the torsion 
tensor is taken, is due to the fact that all other components of this tensor do not 
appear at all in the laplacian of (the otherwise too general) V; in other words, 
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H (g, Q) = |(V) 2 = ±(V) 2 = \A g + Q, with V given by (10 - 11) and V given 
by (18). In the particular case of dimension 2, this is automatically satisfied. In 
the case we actually have assumed, g is Riemannian, the expression (21) is the 
most general invariant laplacian (with zero potential term) acting on functions 
defined on a smooth manifold. This restriction, will allow us to establish a one-to- 
one correspondance between Riemann-Cartan connections of the form (18) with 
(generalized Brownian) diffusion processes. These metric compatible connections 
we shall call RCW geometries (short for Riemann-Cartan- Weyl), since the trace- 
torsion is a Weyl 1-form [10]. Thus, these geometries do not have the historicity 
problem which lead to Einstein's rejection of the first gauge theory ever proposed 
by Weyl. We would like to remark that we first encountered these connections 
on developing a pre-symplectic structure for the derivation of the dynamics of 
relativistic massive spinning systems subjected to exterior gravitational fields 
[57]. 

To obtain the most general form of the RCW laplacian, we only need to 
apply to the trace-torsion one-form the most general decomposition of one-forms 
on a smooth compact manifold. This amounts to give the constitutive equations 
of the particular theory of fluctuations under consideration on the manifold M; 
see [22,26,49]. The answer to this problem, is given by the well known de Rham- 
Kodaira-Hodge theorem, which we present now. 

We consider the Hilbert space of square summable oj of smooth differential 
forms of degree k on M, with respect to vol s . We shall denote this space as L 2 ' fc . 
The inner product is 

«u, (f>»:= / < u(x),<f>(x) > voL, (22) 
Jm 

where the integrand is given by the multiplication between the components 
Wai...a k of w an d the conjugate tensor:g ai/31 . . . g ak ^ k 4>/3 1 ...i3 k ; alternatively, we can 
write in a coordinate independent way: < to(x), 4>{x) > vol g = lo(x) A *<f>(x), with 
* the Hodge star operator, for any oj,<p G L 2,k . 

The de Rham-Kodaira-Hodge operator on L 2,h is defined as 

A fc = -(d + 5) 2 = -(dS + Sd), (23) 

where 5 is the formal adjoint defined on L 2,k+1 of the exterior differential operator 
d defined on L 2 ' k : 

« 5<p,uj >>=<< fadw >>, 
for cp € L 2 ' fc+1 and uj G L 2 ' fc . Then, 5 2 = 0. 
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Let R : (TM ® TM) TM -> TM be the (metric) curvature tensor defined 
by: (V 9 ) 2 Y(v 1 ,v 2 ) = (V 9 ) 2 Y(v 2 , v ± ) + R(v u v 2 )Y(x). From the Weitzenbock 
formula [14] we have 

Ai^(v) = trace (V 9 ) 2 </>(- -) - Ric x (v, 4), 

for v € T X M and Ric x (v\,v 2 ) = trace < R(—,vi)v 2 , — > x . Then, Aq = (V s ) 2 = 
A g so that in the case of k = 0, the de Rham-Kodaira operator coincides with 
the Laplace-Beltrami operator on functions. 

The de Rham-Kodaira-Hodge theorem states that L 2 ' 1 admits the following 
invariant decomposition. Let oj G L 2 ' 1 ; then, 

cj = d f + A 1 + A 2 , (24) 

where / : M — > R is a smooth function on M, Ai is a co-closed smooth 1-form: 
<L4i = —divgAi = 0, and A 2 is a co-closed and closed smooth 1-form: 

5A 2 = 0, dA 2 = 0. (25) 

Otherwise stated, A 2 is an harmonic one-form, i.e. 

AiA 2 = 0. (26) 

Furthermore, this decomposition is orthogonal in L 2,1 , i.e.: 

<< df,A 1 >>=<< df,A 2 >>=<< Ai,A 2 »= 0. (27) 

Remark 1. Note that A\ + A 2 is itself a co-closed one- form. If we con- 
sider an augmented configuration space R x M for an incompressible fluid, this 
last decomposition will be the fluid's velocity. If we consider instead a four- 
dimensional Lorentzian manifold provided with a Dirac-Hestenes spinor operator 
field (DHSOF), one needs the whole decomposition (24) associated to an invari- 
ant density p of the diffusion (i.e. a solution of the equation HQ{g,Q)^p = 0) 
to describe two electromagnetic potentials such that when restricted to the spin- 
plane of the DHSOF, they enforce the equivalence between the Dirac-Hestenes 
equation for the DHSOF on a manifold provided with a RCW connection, and 
the free Maxwell equation on the Lorentzian manifold; see [26, 64]. 

3 Generalized Laplacians on Differential Forms 

Consider the family of zero-th order differential operators acting on smooth 
Worms, i.e. differential forms of degree k (k = 0, . . . , n) defined on M: 

H k (g,Q):=l/2A k + L d , (28) 
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The second term in (28) denotes the Lie-derivative with respect to the vectorfield 
Q. Recall that the Lie-derivative is independant of the metricdor any smooth 
vectorfield X on M 

L x = ixd + dix,, (29) 

where %x is the interior product with respect to X: for arbitrary vectorfields 
X\, . . . , Xk-\ and <p a fc-form defined on M, we have (ix<j>)(Xi, . . . ,X k -±) = 
4>(X, X±, . . . , Xk-i). Then, for / a scalar field, ixf = and 

Lxf = (ixd + di x )f = ixdf = g(X, df) = X(f). (30) 

where X denotes the 1-form associated to a vectorfield X on M via g. We shall 
need later the following identities between operators acting on smooth fc-forms, 
which follow easily from algebraic manipulation of the definitions: 

dA k = A k+1 d, k = 0, . . . , n, (31) 

and 

5A k = A fc _i<5, k = l,...,n, (32) 
and finally, for any vectorfield X on M we have that dLx = Lxd and therefore 
dH k (g,Q) = H k+1 (g,Q)d, k = 0,...,n. (33) 
In (28) we retrieve for scalar fields (k = 0) the operator H(g, Q) defined in 

(21). 

Proposition 1: Assume that g is non-degenerate. There is a one-to-one 
mapping 

V^H k (g,Q) = l/2A fe + L 

between the space of ^-compatible affine connections V with Christoffel coeffi- 
cients of the form 

T% = |j 7 } + ^ {81 Q 7 - 9fh Q a } (34) 

and the space of elliptic second order differential operators on /c-forms (k = 
0, . . . , n) with zero potential term. 
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4 Riemann-Cartan-Weyl Connections and the Lapla- 
cians for Differential Forms 



In this section we shall construct the diffusion processes for scalar fields. 

In the following we shall further assume that Q = Q(r, x) is a time-dependant 
1-form, so that we have a time-dependant RCW connection on M, which we think 
of as a space manifold. The stochastic flow associated to the diffusion generated 
by Ho(g,Q) has for sample paths the continuous curves r i— > x T G M satisfying 
the Ito invariant non-degenerate s.d.e. (stochastic differential equation) 

dx(r) = X{x{t))(1W{t) + Q(t, x(r))dr. (35) 

In this expression, the diffusion tensor X = (X^(x)) is a linear surjection X(x) : 
R m -> T X M satisfying X*X? = g af3 , and {W(t),t > 0} is a standard Wiener 
process on R n . Thus < W T >= and < W*W£ >= 5ijT, where < — > denotes 
expectation with respect to the zero-mean standard Gaussian function on R m 
(to > n). Here r denotes the time-evolution parameter of the diffusion (in a 
relativistic setting it should not be confused with the time variable), and for 
simplicity we shall assume always that r > 0. Consider the canonical Wiener 
space £1 of continuous maps u> : R — > R n , u(0) = 0, with the canonical realization 
of the Wiener process W(t)(uj) = u(t). The (stochastic) flow of the s.d.e. (35) 
is a mapping 

F T : M M, r > 0, (36) 

such that for each oj G the mapping F (. , ui) : [0, oo) x M —>■ M, is continuous 
and such that {F T (x) : r > 0} is a solution of equation (35) with Fq(x) = x, for 
any x G M. 

Let us assume in the following that the components X%, Q a , a, (3 = 1, . . . , n 
of the vectorfields X and Q on M in (35) are predictable functions which further 
belong to C™' € (0 < e < 1, to a non-negative integer), the space of Hoelder 
bounded continuous functions of degree to > 1 and exponent e, and also that 
Q a (r) G L 1 (i?), for any a = l,...,n. With these regularity conditions, if we 
further assume that {x(r) : r > 0} is a semimartingale (see Appendix) on a 
probability space (£l,F,P), then it follows from Kunita [24] that (35) has a 
modification (which with abuse of notation we denote as) 

F t (uj):M^M, F t (lo)(x) = F T (x,u), (37) 

which is a diffeomorphism of class C m , almost surely for r > and w G fi. We 
can obtain an identical result if we assume instead Sobolev regularity conditions. 
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Indeed, assume that the components of a and Q, erf G H S+2 (T* M) and Q@ G 
H S+1 (T*M), 1 < i < m, I < (3 < n, where the Sobolev space H S (T*M) = 
W 2 ' S (T*M) with s > § + m [52]. Then, the flow of (35) for fixed co defines a 
diffeomorphism in H S (M,M) (see [53]), and hence by the Sobolev embedding 
theorem, a diffeomorphism in C m (M, M) (i.e. a mapping from M to M which is 
m-times continuously differentiable as well as its inverse.) In any case, for 1 < m 
we can consider the Jacobian (velocity) flow of {x T : r > 0}. It is a random 
diffusion process on TM, the tangent bundle of M. 

Remarks 2: In the differential geometric approach -pioneered by V. Arnold- 
for integrating NS on a smooth manifold as a perturbation (due to the diffusion 
term we shall present below) of the geodesic flow in the group of volume preserving 
diffeomorphisms of M (as the solution of the Euler equation), it was proved 
that under the above regularity conditions on the initial velocity, the solution 
flow of NS defines a diffeomorphism in M of class C m ; see Ebin and Marsden 
[9]. The difference of this classical approach with the one presented here, is 
to integrate NS through a time-dependant random diffeomorphism associated 
with a RCW connection. As wellknown, these regularity conditions are basic 
in the usual functional analytical treatment of NS pioneered by Leray [45] (see 
also Temam [7]), and they are further related to the multifractal structure of 
turbulence [41]. This diffeomorphism property of random flows is fundamental 
for the construction of their ergodic theory (provided an invariant measure for 
the processes exists), and in particular, of quantum mechanics amd non-linear 
non-equilibrium thermodynamics [10,21,22,26,49]. 

Let us describe now the Jacobian flow. We can describe it as the stochastic 
process on the tangent bundle, TM, given by {v(t) := T xo F T (v(0)) G T Ft ^M, 
v(0) G T XQ M}; here T Z M denotes the tangent space to M at z and T XQ F T is the 
linear derivative of F T at xq. The process {v t ,t > 0} can be described (see [27]) 
as the solution of the invariant Ito s.d.e. on TM: 

dv{r) = V 9 Q(r, v{t))<It + V 9 ' X(v(t))cIW{t) (38) 

If we take U to be an open neighborhood in R n so that TU = U x R n , then 
v(t) = (x(t),v(t)) is described by the system given by integrating (35) and the 
covariant Ito s.d.e. 

dv(T)(x(r)) = V 9 X(x(t)){v(t))(IW(t) + V 9 Q(r, x(r))(u(T))dr, (39) 

with initial condition v(0) = vq € T x ^y Thus, {v(t) = (x(t),v(t)),t > 0} 
defines a random flow on TM. 

Theorem 2 : For any differential 1-form (ft of class C 1,2 (R x M) (i.e. in 
a local coordinate system cj) = a a (T)dx a , with a Q (r, .) G C 2 (M) and a a (.,x) G 
C^i?)) we have the Ito formula (Corollary 3E1 in [27]): 
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4>{v T ) = cP(v )+ j\(V 9 X(v s )dW s + £[£- + H 1 (g,Q)]<l>(v 8 )d8 

+ / T Vmx{x)dW s ){v s ) 
Jo 

+ f T trace d(j)(X(x s )-,V 9 X(v s ))(-)ds (40) 
JO 

In the last term in (40) the trace is taken in the argument — of the bilinear 
form and further we have the mappings 

V 9 Y : TM -> TM; V 9 <p : TM -» T*M. 

Remarks 3 : From (40) we conclude that the infinitesimal generators (i.g., 
for short in the following) of the derived stochastic process is not d T + H\(g, Q), 
due to the last term in (40). This term vanishes identically in the case we shall 
present in the following section, that of gradient diffusions. An alternative method 
which bypasses the velocity process is the construction of the generalized Hessian 
flow further below. Both methods will provide for the setting for the integration 
of the Navier-Stokes equations. 

5 Riemann-Cartan-Weyl Gradient Diffusions 

Suppose that there is an isometric embedding of an n-dimensional com- 
pact orientable manifold M into a Euclidean space R m :f : M — > R m ,f(x) = 
(f 1 (x), . . . , f m (x)). Suppose further that X(x) : R m — ► T X M, is the orthogonal 
projection of R m onto T X M the tangent space at i to M, considered as a subset 
of R m . Then, if e±, . . . , e m denotes the standard basis of R m , we have 

X = X%, with X i = grad f ,i = 1, . . . ,m. (41) 

The second fundamental form [25] is a bilinear symmetric map 

a x : T X M x T X M -> v x M, x € M, (42) 

with v x M = (TxM) 1 - the space of normal vectors at x to M. We then have the 
associated mapping 

A x : T X M x v x M -» T X M, < A x (u, (),v > Rm =< a x (u,v),( > R m, (43) 

for all C £ v x M, u, v € T X M. Let Y(x) be the orthogonal projection onto v x M 

Y(x) = e-X(x)(e),xeM,eeR m . (44) 
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Then: 



V 9 X(v){e) = A x (v,Y(x)e),v G T x M,x G M. (45) 

For any x G M, if we take e\, . . . , e m to be an orthonormal base for R m such that 
ei,. . . ,e m G T X M, then for any v G T X M ,we have 

either V 9 X(u)ej = 0, or X(x)e; = 0. (46) 

We shall consider next the RCW gradient diffusion processes, i.e. for which 
in equation (35) we have specialized taking X = grad/. Let {v T : r > 0} be the 
associated derived velocity process. We shall now give the Ito-Elworthy formula 
for 1-forms. 

Theorem 3. Let / : M — > R m be an isometric embedding. For any differ- 
ential form <j) of degree 1 in C 1,2 (R x M), the Ito formula is 

<f>(v T ) = Hvo)+ [ V 9 <t>{X{x s )dW s )v s + [ T c/>(A x {v s ,Y{x s )dW s ) 
Jo Jo 

+ £[£- + H 1 (g,Q)]<l>(v 9 )ds, (47) 

i.e. d T + Hi(g, Q), is the i.g. (with domain the differential 1-forms belonging to 
C^ 2 {R x M)) of K : r > 0}. 

Proof: It follows immediately from the facts that the last term in the r.h.s. 
of (40) vanishes due to (46), while the second term in the r.h.s. of (40) coincides 
with the third term in (47) due to (45). 

Consider the value & x of a /c-form at x G M as a linear map: $ x : A k T x M — > 
R. In general, if E is a vector space and A : E — > E is a linear map, we have the 
induced maps 

A k A : A k E -> A fe £, A k (v 1 A...A v k ) := Av 1 A ... A 

and 

(dA fe )A : A k E -» A^, 
fe 

(dA fc )i4(u 1 A ... A w fe ) := ^ f 1 A ... A t^'" 1 A ^ A v** 1 A . . . A v k . 

3=1 

For fc = 1, {dA)A = A^4. The Ito formula for fc-forms, 1 < k < n, is due to 
Elworthy (Prop. 4B [27]). 
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Theorem 4. Let M be isometrically embedded in R m . Let Vq G A k T Xo M. 
Set V T = A k (TF T )(V ) Then for any differential form cf) of degree k in C 1 ' 2 ( J RxM), 
1 < jfe < n, 

4>{v T ) = c/>(v ) + [ vmx{x s )dws){v s ) 

Jo 

+ r<A((dA) fc A :Cs (-,y(x s )^)(^)) 

JO 

+ £[^ + H k (g,Q)]<l>(V,)ds (48) 

i.e., <9 r + Hk(g, Q) is the i.g. (with domain of definition the differential forms of 
degree k in C^ 2 {R x M)) of {V T : r > 0}. 

Remarks 4: Therefore, starting from the flow {F T : r > 0} of the s.d.e. 
(35) with i.g. given by d T + Ho(g,Q) , we obtained that the derived velocity 
process {v(t) : r > 0} given by (38) (or (35) and (39)) has H\(g,Q) as i.g.; 
finally, if we consider the diffusion of differential forms of degree k > 1, we get 
that d T + Hk(g,Q) is the i.g. of the process A fc w(r), i.e. the exterior product of 
degree k (k = 1, . . . , n) of the velocity process. In particular, d T + H2(g, Q) is the 
i.g. of the stochastic process {v(t) A v(t) : r > 0}. 

Note that consistently with the notation we have that {A°v T : r > 0} is the 
position process {x T : r > 0} untop of which {A k v T : r > 0}, (1 < k < n) is 
fibered (recall, A°(M) = M). We can resume our results in the following theorem. 

Theorem 5. Assume M is isometrically embedded in R m . There is a one 
to one correspondance between RCW connections V determined by a Rieman- 
nian metric g and trace-torsion Q with the family of gradient diffusion processes 
{A k v T : t > 0} generated by Hk(g, Q), k = 0, . . . , n 

Finally, we are now in a situation for presenting the solution of the Cauchy 
problem 

^ = H k (g,Q T )(x)4>,re[0,T} (49) 
with the given initial condition 

0(0, x) =0 o (x), (50) 

for (j) and <fio fc-forms on a smooth compact orientable manifold isometrically 
embedded in R m . From the Ito-Elworthy formula follows that the formal solution 
of this problem is as follows: Consider the diffusion process on M generated by 
Ho(g,Q): For each r € [0, T] consider the s.d.e. (with s G [0, r]): 

dx T s ' x = X(x T s ' x )dW s + Q{t - s, x T s ' x )ds, (51) 
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with initial condition 

Xq x = x, (52) 
and the derived velocity process {vl~' v ^ = x T s > x , v T s ' v<yX \ < s < r}: 

di r,v(x) = v g X(x T s ' x )(v T s M ^)dW s + V 9 Q(t - s,x T s ' x )(v T s Mx) )ds, (53) 
with initial condition 

vl' v{x) =v{x). (54) 
Then, the C 1 ' 2 (formal) solution of the Cauchy problem defined in [0, T] x M is 
<Mt,*)(AMx)) = ^[^o«' x )(A fc ^^]. (55) 



6 The Navier- Stokes Equation and Riemann-Cartan- 
Weyl Diffusions 

In the sequel, M is a compact orientable (possible with smooth boundary 
dM) n-manifold with a Riemannian metric g. We provide M with a 1-form whose 
de Rham-Kodaira-Hodge (RKH for short) decomposition is 

Q(x) = df(x) + u(x), 5u = —div{u) = 0, 

where / is a scalar field and u is a coclosed 1-form, weakly orthogonal to df, i.e. 
/ g(df, u)vol g = 0. We shall assume that u(x, 0) = u(x) is the initial velocity 1- 
form of an incompressible viscous fluid on M, and that we further have a 1-form 
Q(x,t) = Q a (x,T)dx a whose RKH decomposition is: 

Q(x,t) = df(x,r) + u(x,t), (56) 

with 5u T {x) = 5u(x,t) = (incompressibility condition), and 

J g(dfr,u T )vol(g) = 0, 

which satisfies the evolution equation on M x R (Eulerian representation of the 
fluid): 

^ + V|Q a = -Q^y + fAiQ a , (57) 
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Here v is the kinematical viscosity. In the above notations and in the follow- 
ing, all covariant operators act in the M variables only. In the formulation of 
Fluid Mechanics in Euclidean domains, Q(x,t) receives the name of (Buttke) 
"magnetization variable" [1]. 

Remarks 5: We recall that to take the RKH decomposition of the ve- 
locity of a viscous fluid is a basic procedure in Fluid Mechanics [1,6,7,9]. We 
shall see below that Q and in particular u are related to a natural RCW geom- 
etry of the incompressible fluid. In the formulation of Quantum Mechanics and 
of non-linear non-equilibrium thermodynamics stemming from RCW diffusions, 
we have a RKH decomposition of the trace-torsion associated to a stationary 
state; see [10,21,22,26,49]. This decomposition allows to associate with the diver- 
genceless term of the trace-torsion a probability current which characterizes the 
time-invariance symmetry breaking of the diffusion process, and is central to the 
construction of the ergodic theory of these flows. 

Equation (57) is the gauge-invariant form of the NS for the velocity 1-form 
u(x,t). Indeed, if we substitute the decomposition of Q(x,t) = Q T (x) into (57) 
we obtain, 

£ + Viu T = vA lUr - d(-L + V| r / + -K| 2 - uA g f). (58) 

Consider the operator P of projection of 1-forms into co-closed 1-forms: Puj = a 
for any one-form oj whose RKH decomposition is uj = df + a, with 5a = 0. From 
(32) we get that 

PA lUr = Antr, (59) 

and further applying P to (57) we finally get the well known covariant NS (with 
no exterior forces; the gradient of the pressure term disappears by projecting with 
P [1,9]) 

du 

— + P[VIu t ]-uA iUt = 0. (60) 

Conversely, starting with equation (58) which is equivalent to NS we obtain (57). 
Note that Q T and u T differ by a differential of a function for all times. Multiplica- 
tion of (58) by / — P (I the identity operator) yields an equation for the evolution 
of / which is only arbitrary for r = 0. Now we note that the non-linearity of NS 
originates from applying P to the term 

V| T «r = iu T du T , 

which taking in account (29) can still be written as 

Lu T u T - diu T u T = L^Ur - l/2d(\u T \ 2 ). (61) 
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Applying P to (61), we see that the kinetic energy term there disappears and the 
non-linear term in NS can be written as 



P[Vlu T ]=P[h ilT u T ]. (62) 
Therefore, from (28) and (62), NS (60) takes the final concise form 

^=PH 1 (2vg,^- UT )u T . (63) 
or 2v 

Therefore we have found that NS for the velocity of an incompressible fluid is a 
a non-linear diffusion equation. 

Let us introduce the vorticity two-form 

Q T = du T . (64) 

Note that also £l T = dQ T . Now, if we know £l T for any r > 0, we can obtain u T (or 
still Q T ) by inverting the definition (64). Namely, applying 5 to (64) and taking 
in account (23) we obtain the Poisson-de Rham equation (would g be hyperbolic, 
it is the Maxwell-de Rham equation [10a]) 

Aiu T = -Sn T . (65) 

and an identical equation for Q T . (Note that if we know Q T we can reconstruct 
f T by solving — 5Q T = div(Q T ) = A g f T , for any r.) From the Weitzenbock 
formula, we can write (65) showing the coupling of the Ricci metric curvature to 
the velocity u = u a (x,T)dx a : 

(V 9 ) V - R a p U l?dx a = -Sn T . (66) 

with R^,{g) = Raa {9)1 the Ricci (metric) curvature tensor. Thus, the vorticity 
S1 T is a source for the velocity one-form u T , for all r; in the case that M is a 
compact euclidean domain, equation (66) is integrated to give the Biot-Savart 
law of Fluid Mechanics [1,39]. 

Now, apply d to (63) and further RKH decompose L_^ T u T = a T + dp T (with 
p T the pressure at time r); in account that 

dPL-^Ur = da T = d(a T + dp T ) = dL^Ur = L_,~ aT du T = L_„ T T , 

and that from (31) we have that dA\u T = A2^ T; we therefore obtain the linear 
evolution equation 

= H 2 {2ug, —u T )n T . (67) 
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Thus, we have proved that the Navier-Stokes equation is a linear diffusion 
equation generated by a RCW connection. This connection has 2ug for the 
metric, and the time-dependant trace-torsion of this connection is —uj2v. Then, 
the drift of this process does not depend explicitly on u, as it coincides with the 
vectorfield associated via g to — u T , i.e.—u T . Notice that when v tends to zero, 
i.e. the Euler equation, the trace-torsion becomes singular. 

Theorem 6 : Given a compact orientable Riemannian manifold with metric 
g, the Navier-Stokes equation (63) for an incompressible fluid with velocity one- 
form u = u(t,x) such that Su T = 0, assuming sufficiently regular conditions, 
are equivalent to a linear diffusion process for the vorticity given by (67) with 
u T satisfying the Poisson-de Rham equation (65). The RCW connection on M 
generating this process is determined by the metric 2vg and a trace-torsion 1-form 
given by —uj2v. 

Remarks 6: We would like to recall that in the gauge theory of gravitation 
[46,57] the torsion is related to the translational degrees of freedom present in 
the Poincare group, i.e. to the gauging of momentum. Here we find a similar, 
yet dynamical situation, in which the trace-torsion is related to the velocity. 
We would like to point out further, than on setting the skew-symmetric torsion 
to be zero on taking RCW connnections, due to the fact that only the trace- 
torsion appears in the laplacian generating generalized Brownian motions, we are 
setting to zero the inertial fields which can be associated with the skew-symmetric 
torsion. Thus, generalized Brownian motions are not generated by inertial fields 
[64]. 

7 Random Diffeomorphims and the Navier-Stokes Equa- 
tions 

In the following we assume additional conditions on M, namely that / : 
M — > R m is an isometric embedding, and that M has no boundary. 

Let u denote a solution of (63) (which we assume that exists for r G [0,T]), 
and consider the flow {F" : r > 0}) of the s.d.e. whose i.g. is 
from (35) we know that this is the flow defined by integrating the non-autonomous 
Ito s.d.e. with X = V/) 

dx"' T ' x = [2u]^X(x u ' T ' x )dW(T) - u(r,x u ' T ' x )dT,x u '°' x = x,t G [0,T]. (68) 

We shall assume in the following that the diffusion tensor X and the drift 
u T have the regularity conditions stated in Section 4, so that the randoms flows 
of (68) is a diffeomorphism of M of class C m , m > 1. 
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Theorem 7: Equation (68) is a random Lagrangian representation for the 
fluid particles positions, i.e. x(r) is the random position of the fluid particles of 
the incompressible fluid whose velocity obeys (63). 

Remark 7: Note that the drift of the derived process {v T : r > 0} is 
minus the deformation tensor of the fluid (see equation (71) further below). This 
will have a crucial role in the solution for the vorticity equation as well as the 
kinematic dynamo problem. We further note that if in (68) we set the viscosity 
to zero, we get the classical flow of the Euler equation. 

7.1 Cauchy Problem for the Vorticity 

Let us solve the Cauchy problem for 0(r, x) of class C m in [0,T] x M 
satisfying (67) with initial condition Qq(x). 

For each r G [0, T] consider the s.d.e. (with s G [0, r]) (obtained by running 
backwards the Lagrangian representation (68) above): 

dx u s ' T ' x = (2u) 1 2X(x" s ' T ' x )dW s - u{t - s,x" s ' T ' x )ds, (69) 
with initial condition 

xT X = *, (70) 
and the derived velocity process {v U s ,T ' v ^ = (Xg' T ' x , Vs' T ' v< " x ^), < s < r}: 
dv u s ' T ' v{x) = {2u)\v 9 X{x v s ^ x ){v v s ^ v ^)dW s - V 9 u{t - s,x v s ' T ' x )(v»' T ' v W)ds, (71) 
with initial condition 

v% T ' v W = v(x). (72) 

Theorem 8: Let Q T (x) be a bounded C 1,2 solution of the Cauchy problem; 
then it follows from the Ito formula-Elworthy (48) (with k = 2) is 

Q T (A 2 v(x)) = E x [n (x^> x )(A 2 v»> T > v W)} (73) 

where the expectation value at x is taken with respect to the measure on the 
process {x" ,T,x : r G [0,T]} (whenever it exists): 

Proof: It follows just from applying the Ito-Elworthy formula for 2-forms.c.q.d. 

Remarks 8: We would like to examine the physical interpretation of the 
representation (73). To determine the vorticity at time r on a point x evaluated 
on a bivector A 2 v(x), we run backwards the random lagrangian representations 
starting at time r at x, and its Jacobian flow starting at v(x) along which we 
transport the time-zero vorticity, and further we take the mean value along all 
possible random paths. Furthermore, the transport of f^o along this Jacobian 
flow, indicates that the original vorticity is acted upon by the fluid-deformation 
tensor and the gradient noise term. 
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8 Integration of the Poisson-de Rham equation 



In (68) we have that u T verifies (65), for every r > which from (28) we can 
rewrite as 

Hi(g,0)u T = --dn T , for any r > 0. (74) 

This last representation together with the Ito-Elworthy formula for 1-forms, in- 
dicates automatically what the representation for the solution is: We have to 
construct a Jacobian process on TM originated from the derivative of the scalar 
diffusion with zero drift and diffusion tensor defined by g. 

Thus, consider the autonomous s.d.e. generated by Ho(g,0) = ^A 9 : 

dx 9 ' x = X(x 9 s )dW s ,x 9 ' x = x. (75) 

We shall solve the Dirichlet problem in an open set U (of a partition of unity) 
of M given by (74) with the boundary condition u T = <fi on dU, with <p a given 
1-form. Then one can "glue" the solutions and use the strong Markov property 
to obtain a global solution (cf. [31]). Consider the derived velocity process 
W 9( s ) = (x 9 (s),v 9 (s)) on TM, with v 9 (s) G T x3{s) M, whose i.g. is Hi(g,0), i.e. 
from (35) we have: 

dvf v{ - x \xf x ) = V 9 X(x 9 ' x (s)){v 9 ' v ^(s))dW(s), (76) 

with initial velocity v 9 ' v ( x \0) = v(x). Notice that equations (75&76) are obtained 
by taking u = in equations (69&71) respectively, and further rescaling by 
(2z/)~2 . Then if u T is a solution of (74) for any fixed r € [0, T], applying to it the 
Ito-Elworthy formula and assuming further that 5Q T is bounded, we then obtain 
that the formal C 1,2 solution of the Dirichlet problem is given by: 

Ur(x){v(x)) = E%[<j>[&?)W?W + fj hn T (xf X (vr^)d S ] 

= J \4>{x 9 ?){v 9 ^ e) + l/2 J T J sn T (y)(vl^(y)ds]p 9 (s,x,y)vol g (y), (77) 

where r e = inf{s : x 9,x ^ U}, the first-exit time of U of the process {x 9 ' x }, 
and E B denotes the expectation value with respect to p 9 (s,x,y), the transition 
density of the s.d.e. (76), i.e.the fundamental solution of the heat equation on 
M: 

d T p(v) = H (g, 0)(y)p(y) = l/2A g p(y) (78) 
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with p(s, x, — ) = S x as s J. 0. 

Theorem 9: Assume 5 is uniformly elliptic, and U has a C 2,e -boundary, 
and furthermore g a @ and 50 T are Hoelder-continuous of order e on J7 and u T 
is uniformly Hoelder-continuous of order e, for r € [0,T]. Then the solution of 
the Dirichlet problem above has a unique solution belonging to C 2 ' e (U) for each 
r £ [0,T] (the maximum principle) [31,47]. Assume instead that u T £ ff 1 (T*L r ) 
for each r £ [0,T], i.e. belongs to the Sobolev space of order 1. If 5Q T E 
i7 fe_1 (A 1 (T* [/")), then u T E F fc+1 (A 1 (T*C/)), for fc > 1 and r € [0, T] (cf. [53]). 

Notice that in the representations (73&77), the local dependance on the 
curvature is built-in (the curvature is defined by second-order derivatives). This 
dependance might be exhibited through the scalar curvature term in the Onsager- 
Machlup lagrangian appearing in the path-integral representation of the funda- 
mental solution of the transition densities of equations (73) and (77) [35,44]. 
There is further a dependance of the solution on the global geometry and topol- 
ogy of M appearing through the Riemannian spectral invariants of M in the 
short-time asymptotics of these transition densities [28,29,43]. 

9 Kinematic Dynamo Problem of Magnetohydrody- 
namics 

The kinematic dynamo equation for a passive magnetic field transported by 
an incompressible fluid, is the system of equations [56] for the time-dependant 
magnetic vectorfield B(t,x) = B T (x) on M defined by ib t ^{x) = oj t (x) (for 
t > 0), satisfying 

dfio + (L{i T - z/^A^Vt = 0, w(0, x) = u(x), < t, (79) 

where v m is the magnetic diffusivity, and we recall that (i = vol(g) = det(g) ^dx l A 
...Adx n is the Riemannian volume density ((x 1 , . . . , x n ) a local coordinate system 
on M), and uj G A™~ 1 (i? x T*M). In (79), u is assumed given, and it may either 
be a solution of NS, or of the Euler equation given by setting v = in (63). From 
the definition follows that divB T = 0, for any r > 0. Now we note that from (28) 
we can rewrite this problem as 

8 T u; = H n ^(2u m g,-^u T )u; T ,u;(0,x) =lo(x),0<t, (80) 

as a linear evolution equation for a (n — l)-form, similar to the evolution Navier- 
Stokes equation for the vorticity. Now if we assume that there is an isometric 
embedding / : M — > R m , so that the diffusion tensor X = V/, we can take the La- 
grangian representation for the scalar diffusion generated by Ho(2i> m g, —-^p^u T ), 
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i.e. the Ito s.d.e. given by substituting v m instead of v in equation (68), and we 
consider the jacobian process given by (69 — 72) and with u m instead of u, then 
the formal C 1 ' 2 solution of (80) defined on [0, T] x M for some T > 0, is given by 

Co T (K n - l v{x)) = E x [u (x'T' T > x )(A n - 1 v» m ' T > v M)} (81) 

Remarks 8: We note that similarly to the representation for the vortic- 
ity, instead of the initial vorticity, now it is the initial magnetic form which is 
transported backwards along the scalar diffusion, where now the parameter is the 
magnetic diffusivity, and along its way it is deformed by the fluid-deformation 
tensor and the gradient diffusion tensor noise term (this accurately represents the 
actual macroscopical physical phenomena yet in a microscopic approach) , and fi- 
nally we take the average for all those paths starting at x. For both equations 
as well as the Poisson-de Rham equation, we have a microspic description which 
clearly evokes the Feynman approach to Quantum Mechanics through a sum- 
mation of the classical action of the mechanical system along non-differentiable 
paths. In distinction with the usual Feynman approach, these Brownian integrals 
are well defined and they additionally have a clear physical interpretation which 
coincides with actual experience. 

10 Random Implicit Integration Of The Navier-Stokes 
Equations For Compact Manifolds 

Up to this point, all our constructions have stemmed from the fact that for 
gradient diffusion processes, the Ito-Elworthy formula shows that the random 
process on A 2 TM given by {A 2 v T : r > 0} with {v T : r > 0} the jacobian 
process fibered on the diffusion process {A°v T = x T : r > 0} on M given by 
(68) , is a random Lagrangian flow for the Navier-Stokes equation. Our previous 
constructions have depended on the form of the isometric embedding of M. This 
construction is very general, since that from a well known theorem due to J. Nash 
(1951), such an immersion exists of class C 1 for any smooth manifold (cf. [53]). 
(Furthermore, our assumption of compactness is for the obtention of a random 
flow which is defined for all times, and gives a global diffeomorphism of M. The 
removal of this condition, requires to consider the random flow up to its explosion 
time, so that in this case we have a local diffeomorphism of M.) 

There is an alternative construction of diffusions of differential forms which 
does not depend on the embedding of M in Euclidean space, being thus the 
objective of the following section its presentation. A fortiori, we shall apply these 
constructions to integrate NS and the kinematic dynamo problem. 
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10.1 The Generalized Hessian Flow 



In the following M is a complete compact orientable smooth manifold with- 
out boundary. We shall construct another flow in distinction of the derived 
flow of the previous sections, which depends explicitly of the curvature of the 
manifold, and also of the drift of the diffusion of scalar s. We start by consid- 
ering an autonomous drift vector field Q (further below we shall lift this condi- 
tion) and we define a flow W k > Q on A k T*M (1 < k < n) over the flow of (35), 
{F t (xq) : t > 0}, by the covariant equation 

a T (V ) = -l/2R k (W k ^(V )) + (dA k )(V 9 Q(.))(W k ^(V )), (82) 

OT 

with V £ A k T XQ M; in this expression the operator 4^ denotes parallel transport 
along the curves x T with V 9 ; R k is the Weitzenbock term (see [14]) appearing in 
the Weitzenbock formula for fc-forms: = (V s ) 2 — R k . Let be a fc-form in 
C 1,2 (M), and V T = W k ' Q (V ) and x T = F T (x ); then from the Weitzenbock and 
general Ito formula we have the following Ito-Elworthy formula [27]: 

4>{V T ) = <t>{V ) + [ T VmX(x s )dW s )V s + [ T [H k (g, QMV s )ds; (83) 
Jo Jo 

In other words, H/,(g, Q) is the i.g. of {V T }. In the case that Q is exact the flow 
V T is called the Hessian flow. Assume that l/2R k - (dA) k (V g Q))(.)) is bounded 
below, i.e. for any V € A k TM with |V| = 1 we have 

-oo< C k (Q)=inf^ < R k {V),V > - < (dA) k (V 9 Q)(.))V, V >, (84) 

where we have denoted by <, > the induced metric on A k TM. 

Proposition 2 (Elworthy [27]) Assume that \R k - (dA) k (V 9 Q)(.)) is 
bounded below. Define P k : L°°A k T*M -► L°°A k T*M by 

P k ( ( f>)(V) = E( ( f>(W k ^(V)) (85) 

for V G A k T x M, cf) G L°° A k T* M . Then {P k : r > 0} is a contraction semigroup 
of bounded continuous forms and is strongly continuous there with i.g. agreeing 
with H k (g,Q) on C 2 (M). 

Under the above conditions we can integrate the heat equation for bounded 
twice differentiable fc-forms of class C 2 (0 < k < n) and in the general case of a 
non-autonomous drift vector field Q = Q T (x). Indeed, for every r > consider 
the flow Vg = W k, P over the flow of {x T s : < s < r}, given by the equation 

dxf T = X(xf T )dW s + Qr-s(x x s ' T )ds, x x / = X, 
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obtained by integration of the equation 
D^V 



, h y v ) = -\/2R«{VJ{v )) + (dA«)(V<3Q T - s (.))(V;(v )), (86) 

with vo = Vq G T x M. Then, applying the Ito-Elworthy formula we prove as 
before that if a T is a bounded C 1,2 solution of the Cauchy problem for the heat 
equation for k forms: 

^a T = H k (g,Q)a T (87) 

with initial condition cxq(x) = a(x) a given /c-form of class C 2 , then the solution 
of the heat equation is 

a T (v(x)) = E x [a (V T T (v(x))], 



with V^{x) the generalized Hessian flow over the flow {F T (x) : r > 0} of {x T T : 
t > 0} with initial condition v(x). 

To integrate the Poisson-de Rham equation we shall need to consider the so- 
called Ricci-flow = W^ ,0 (lo) : TM — ► TM over the random flow generated 
by Ho(g, 0), obtained by integrating the covariant equation (so we fix the drift to 
zero and further take k = 1 in (82)) 

—g^M = -^Ric(W^(v ), -),vo G T X0 M (89) 

where Ric : TM TM -> is the Ricci curvature and B«c(u, -) G T^M is the 
conjugate vector field defined by < Ric(v, —),w >= Ric(v,w), w G T^M. 

10.2 Integration of the Cauchy problem for the Vorticity on 
Compact Manifolds 

Theorem 10: The integration of the equation (67) with initial condition 
^(0, ) = Q yields 

n T (v(x)) = E x [n (V T T (v(x))} (90) 

where {V^ : r > 0} is the solution flow over the flow of {x" ,T,x : r > 0} (see 
equation (69)) of the covariant equation 

— ^ (v(x)) = -vR 2 (W*>- U0 (v(x))) 

OT 

- (dA 2 )(V9u (.))(W 2 ^ (v(x))) (91) 

with initial condition v(x) G T X M. In this expression, V 9 uq(.) is a linear trans- 
formation, A, between T*M and T X M, and dA 2 (A) : T X M t\T x M -> T X M t\T x M 
is given by dA 2 A(wi A U2) = A i>2 + v\ A At>2, for any v±, V2 G T X M, x G M. 
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10.3 Integration of the Kinematic Dynamo for Compact Mani- 
folds 

Substituting the magnetic diffusivity v m instead of the kinematic viscosity in (69) 
and we further consider {Vg : s € [0, r]} given by the solution flow over the flow 
of {x v s m,T,x : s G [0,t]} (see equation (69)) of the covariant equation 

-of-{v{x)) = -v™R n -\Vj{v{x))) + (dA^iVSQrsmVJWxm (92) 

with v(x) = V T £ T X M. Then, the formal C 1 ' 2 solution of (80) is 

u T (v(x))=E x [u (V;(v(x))}. (93) 

with E x denoting the expectation valued with respect to the measure on {x^ m,T,x } 
(whenever it exists). 

10.4 Integration of the Poisson-de Rham Equation for the Ve- 
locity 

With the same notations as in the case of isometrically embedded manifolds, we 
have a martingale problem with a bounded solution given by 

u T (v(x)) = E^[4>(W^(v(x))) + 1/2 f e 5n T (W?(v(x)))ds} (94) 

Jo 

11 Representations Of NS On Euclidean Space 



In the case that M is euclidean space, the solution of NS is easily obtained 
from the solution in the general case. In this case the isometric embedding / 
of M is realized by the identity mapping, i.e. f(x) = x, Vx € M. Hence the 
diffusion tensor X = I, so that the metric g is also the identity. For this case 
we shall assume that the velocity vanishes at infinity, i.e. u% — > as \x\ — > oo. 
(This allows us to carry out the application of the general solution, in spite of the 
non-compacity of space). Furthermore, r e = oo. The solution for the vorticity 
equation results as follows. We have the s.d.e. (see (69) where we omit the 
kinematical viscosity, for simplicity) 

dx T s ' x = -u{t - s,x T s ' x )ds + {2v)^dW s ,x T ' x = x,se [0,r]. (95) 

The derived process is given by the solution of the o.d.e. (since in (71) we have 
VX = 0) 

d yT,x,v(x) = _ Vn (r - s, x T s ' x )(v T s ' x ' v W)ds, vl' v{x) = v(x) £R a ,se [0, r], (96) 
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Now for n = 3 we have that the vorticity f2(r, x) is a 2-form on R 3 , or still 
by duality has an adjoint 1-form, or still a function, which with abuse of notation 
we still write as 0(r, .) : -R 3 — > -R 3 , which from (73) we can write as 

n(T,x) = E x [v T T > x > I n (x T T > x )}, (97) 

where E x denotes the expectation value with respect to the measure (if it exists) 
on {x T T ,x : t > 0}, for all x € i? 3 , and in the r.h.s. of (97) we have matrix 
multiplication Thus, in this case, we have that the deformation tensor acts on 
the initial vorticity along the random paths. This action is the one that for 
3D might produce the singularity of the solution. Note that in Euclidean space 
there is no gradient-noise contribution to the folding of the initial vorticity by 
the action of the fluid deformation tensor. 

In the case of R 2 , the vorticity can be thought as a pseudoscalar, since 
fl T (x) = Qri^dx 1 A dx 2 , with £l T : R 2 — > R, and being the curvature identically 
equal to zero, the vorticity equation is (a scalar diffusion equation) 

-^ = H (2uI,—u T )n T (98) 
so that for = Q given, the solution of the initial value problem is 

n( T ,x) = E x [n(x T T ' x )} (99) 

This solution is qualitatively different from the previous case. Due to a geometri- 
cal duality argument, for 2D we have factored out completely the derived process 
in which the action of the deformation tensor on the initial vorticity is present. 
Furthermore, the solution of equation (75) is (recall that X = I) 

x 9f = x + w r , (100) 

and since VX = 0, the derived process (see (76)) is constant 

v g,x,v(x) =v ( a; ) j Vr G [0,T]. (101) 

so that its influence on the velocity of the fluid can be factored out in the repre- 
sentation (77). Indeed, we have 

u T (x)(y(x)) = E%[ -Sn T (x + W s )(v^^)(s))ds] 
JO 2 

= E^[j o -Sn T (x + W s )ds(v(x))} 
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for any tangent vector v(x) at x, and in particular (we take v(x) = I) we obtain 



poo 1 



u T (x)=E*[ -5n T (x + W s )ds]. (102) 
jo 2 

In this expression we know from (78) that the expectation value is taken with 
respect to the standard Gaussian function, p g (s,x,y) = (4ns)~exp(— ^ x ^ ). 

Let us describe in further detail this solution separately for each dimension. 
We note first that if Q T £ L 1 n C\ (where C\ means continuously differentiable, 
bounded with bounded derivatives) 

E[5Q T (x + W s )] = 5E[n T (x + W a )\ (103) 

In case n = 2, for a 2-form j3 on M we have 5(3 = 5(f3dx 1 A dx 2 ) = —(dzfidx 1 — 
di(3dx 2 ) = — V -1 -/?. In case n = 3, for a vorticity described by the 1-form (or a 
vector- valued function) Q T : R 3 ^ R 3 adjoint to the vorticity 2-form Q T , we have 
that 

5Q T = -dn T = -rotfV (104) 

Therefore, we have the following expressions for the velocity: When n = 2 we 
have 

roo 1 

u T (x)= --V ± E^[n T (x + W s )]ds (105) 
Jo 2 

while for n = 3 we have 

u T (x) = r —dE^[n T (x + W s )]ds. (106) 
Jo 2 

Now we can obtain an expression for the velocity which has no derivatives of 
the vorticity; this follows the basic idea in a non-geometrical construction given 
by Busnello who starts with the stream function (of an unbounded incompressible 
fluid) instead of the velocity [54]. Consider the semigroup generated by Hq(I, 0) = 
^A, i.e. P s Q T (x) = E[Q T (x + W s )] (in the case n = 3 this means the semigroup 
given on each component of Q). From the Elworthy-Bismut formula valid for 
scalar fields (see [52]) we have that (in the following e«,i = 1,2,3 denotes the 
canonical base in R? or R?) 

diP s & T {x) =< dP s n(x), ei >= -E^[Cl{x + W a ) T < ei ,dW r >] 

s Jo 

= -E*[£l(x + W s ) f dW l r ] = -Ep T (x + W s )W l s }. (107) 
s Jos 



30 



Therefore, for n = 2 we have from (105, 107) 



(*) = "/ 
Jo 



-E^[n T (x + W s )W s ± ]ds 



(108) 



o 2s 



where = (W^W 2 ) 1 = (W s 2 ,— W}). Instead, for n = 3 we have from 
(106, 107) that 



where x denotes the vector product and W = (W 1 , W 2 , W 3 ) G R 3 . 
Thus, we have obtained the representations for NS in 2D and 3D: 

11.1 Integration of the Kinematic Dynamo Problem in Euclidean 



With the notations in this section, the kinematic dynamo problem in 3-D can be 
solved as follows. As for the vorticity, the magnetic field is for n = 3 is a 2-form 
on R 3 , or still by duality has an adjoint 1-form (so the argument turns to work 
out as well for 2D), or still a function, which with abuse of notation we still write 
as lj(t, .) : R 3 — > R 3 , which from (81) we can write as 



where E x denotes the expectation value with respect to the measure (if it exists) 
on {x^ ,x : r > 0}, for all x € R 3 , and in the r.h.s. of (110) we have matrix mul- 
tiplication Thus, in this case, we have that the deformation tensor acts on the 
initial vorticity along the random paths. This action is the one that for 3D pro- 
duces the complicated topology of transported magnetic fields. This solution was 
obtained independently by Molchanov et al [50] and further applied in numerical 
simulations (see Ghill and Childress [51] and references therein). 

12 The Navier-Stokes Equation is Purely Diffusive 
For Any Dimension Other than 1 

12.1 Motivations 

We have given up to now a derivation of diffusion processes starting from 
gauge theoretical structures, and applied this to give implicit representations 
for the covariant Navier-Stokes equations. These constructions were possible as 




(109) 



space 



u(t,x) = E x [v t t > x > i u;o(x t t > x )], 



(110) 
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they stemmed from the extremely tight relation existing between the metric- 
compatible Riemann-Cartan-Weyl connections, and the diffusion processes for 
differential forms, built untop of the diffusions for scalar fields. As we saw already 
this stemmed from the fact that there is a one-to-one correspondance between said 
RCW connections and the scalar diffusion processes {x T : r > 0} with drift given 
by Q and diffusion tensor X. As we can easily check from (34), this construction 
is valid for n ^ 1. This leads to conjecture that in a gauge theoretical setting 
and further applying stochastic analysis, one could do away with the drift, in 
any dimension other than 1. If this would be the case, then we could apply this 
construction to the Navier- Stokes equation, which thus in any dimension other 
than 1 would turn to be representable by random lagrangian paths which do not 
depend explicitly on the velocity of the fluid, since they would be purely diffusive 
processes. 

12.2 More On Connections 

Consider the map M x R m — > TM — > which we assume that it has a right 
inverse Y : TM — > M x R m . Here, Y = is the adjoint of X with respect to 
the Riemannian metric on TM induced by X, Y = X* . Write X{x) = X(x, .) : 
R m -► TM. For u G TM, let Z u G F(TM) defined by 

Z u (x) = X(x)Y(tt(u))u. (Ill) 

Proposition 3: There is a unique linear connection V on TM such that for 
all uq G T XQ M, x G M, we have that 

V VO Z UO = 0. (112) 

It is the pushforward connection defined as 

V V0 Z := X(x )d(Y(.)Z(.))(v ),v G T Xo M, Z G T(TM), (113) 

where d is the usual derivative of the map Y{.)Z{.) : M — ► R m . 

Proof: The above definition defines a connection. Let V be any linear 
connection on TM. We have 

Z{.) = X(.)Y{.)Z{.). (114) 

Then, for v G T XQ M, 
V V Z = X(x )d(Y(.)Z(.))(v) + V v [X(.)(Y(x )Z(x ))} = V V Z + V V Z Z ^\ (115) 
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Since V is a connection by assumption, and since the map 

TM x TM — > TM, (v, u) i — ► V„Z (116) 

gives a smooth section of the bundle Bil(TM x TM;TM), then V is a smooth 
connection on TM. Taking V = V we obtain a connection with the desired 
property. 

Theorem 11: Let Y be the adjoint of X with respect to the induced metric 
on TM by X. Then, V is metric compatible, where the metric is the one induced 
by X on TM, which we denote by g. Moreover, since M is finite-dimensional, 
any metric-compatible connection for any metric on TM can be obtained this 
way from such X and R m . 

Proof: We have 

2g(V v U, U) = 2(X(x )(d(Y(.)U(.))(v), U(x )) 

= 2~g(d(Y(.)U(.))(v), Y(x )U(x )) = d(g(U, U)(v), (117) 

so that V is indeed metric compatible. By the Narasimhan-Ramanan theorem 
on universal connections [51], any metric compatible connection arises likes this. 
Indeed, V is the pull-back of the universal connection over the Grassmanian 
G(m,n) of re-planes in R m by the map x i— » [Imagey(x) : T X M — > R m \. In 
particular, the RCW connections arise from such a construction, c.q.d. 

Two connections, V a and V b on TM give rise to a bilinear map D ab : TM x 
TM -> TM such that 

V^U = V b v U + D ab (V, U), U, V € r(TM). (118) 

Choose V b = V 9 , the Levi-Civita connection of a certain Riemannian metric g. 
Consider 

V V U = V 9 V U + D(V,U), (119) 

where we decompose D into 

D(u,v) =A(u,v) +S(u,v), (120) 

where 

A(u, v) = -A(v, u), S{u, v) = S{v, u). (121) 

Since the torsion tensors T a and T b of any two connections V a and V fe respec- 
tively, are connected through the expression 

T a (u, v) + T b (u, v) = D ab (u, v) - D ab (v, u). (122) 
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which for the case of V b = V 9 as T b = 0, we can write for D ab = D, the identity 
f(u, v) = D(u, v) - D(v, u) + 2A(u, v). (123) 

Thus, 

A(u, v) = - f(u, v), u, v, G T(TM),f = T a . (124) 
Therefore, the decomposition in (119) is written in the form 

V V U = V 9 V U + -f (u, v) + S(u, v),u,v€ r(TM). (125) 

which is nothing else than the original decomposition for a metric compatible 
Cartan connection given in (10&11). As we know already from the decomposition 
(11), we have the following lemma. 

Lemma A connection V is metric-compatible if and only if the map D(v, .) : 
TM — > TM is skew-symmetric for each v G TM, i.e. 

g(D(v, Ul ),u 2 ) + g{D{v, u 2 ), Ul ) = 0, Ul ,u 2 £ T(TM). (126) 

Equivalently, 

g(S(u 1 ,u 2 ),v) = -g(f(v, ui), u 2 ) + -g(f(v, u 2 ), u 2 ). (127) 
Consequently, for Ui, U 2 ,V in T(TM), we have 

g(D(V, U 2 ) = \g{f{V, UJ, U 2 ) + \g{f{U 2 , V), U{) + \g{f{U 2 , U x ), V), 

(128) 

which is decomposition (10). 

Proof: Take V,U U U 2 G T(TM). Then, 

%(C/i, U 2 )){V) = g{V 9 v U u U 2 ) + g{U u V g v U 2 ) 

= g{V 9 v U u U 2 ) + g(U lt V g v U 2 ) - g(D(V, Ui), U 2 ) - g^DiV, U 2 ). (129) 
So, V is metric compatible if and only if 

g(D(V,U 1 ),U 2 )+g(U 1 ,D(V,U 2 )) = 0. (130) 

Now, writing D = A + S we get 

g(A(V, f/i), U 2 ) + g(A(V, U 2 ), U{) = -g(S(V, U{), U 2 ) - g(S(V, U 2 ), U{). (131) 
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We now observe that for an alternating bilinear map L : TM x TM — ► TM, 

Cyl[g(L(v, m), u 2 ) + g{L(v, u 2 ),u l )} = 0, (132) 

where Cyl denotes cyclic sum. Taking the cyclic sum in equation (131) and apply 
(132) to A, we thus obtain Cylg(S(V, Ui), U 2 ) = which on further substituting 
in (131) we obtain 

g(A(V, U x ), U 2 ) + g(A(V, U 2 ), U{) = g(S(U u U 2 ), V). (133) 

12.3 The Trace- Torsion Is Dynamically Redundant in Any Di- 
mension Other Than 1 

Let us return to our original setting of Section 1. We assume a metric-compatible 
Cartan connection, which we now write as V with torsion tensor T. The following 
result is a reduction of a more general result due to Elworthy, Le Jan and Li 61]. 

Theorem 12: Assume M has dimension bigger than 1. Consider the lapla- 
cian on 0- forms Ho(g, Q) where Q is the trace-torsion 1-form of V, 

Q(u) = trace g(f(-, u), -). (134) 

Assume further that we can write the laplacian H p (g, Q) on p-forms (0 < p < n) 
in the Hormander form: 

m 

-^L Vi L Vi + Lz (135) 
z i=i 

where Z is a vectorfield on M, V : M x R m — > TM is a smooth surjection, linear 
in the second variable, and V{ is defined by 

m 

V(x,e) = V(x)e = Y / y i ( x ) < e,ei>, (136) 
i=i 

whith e±,...,e m the standard orthonormal basis for R m . (Since V 9 is metric 
compatible, from Theorem 11 and the transformation rules between Stratonovich 
and Ito calculi, we can always introduce a defining map V for V 9 that gives such 
decomposition (c.f. [27])). Then, there exists a map K : M x R m — > TM linear 
on the second variable, such that the solution to the Stratonovich equation 

dx T = K(x T )odW T , (137) 

has Ho(g,Q) for infinitesimal generator, i.e. 

-. m 

H (g,Q) = -J2L Kt L Kt . (138) 
z i=i 
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In other words, the Ito s.d.e. given by (35) admits a driftless representation given 
by (137). 

Proof: Set for the original drift vectorfield Q (the g-conjugate of Q), the 
decomposition 

z 1=1 

A connection V suitable for this is such that 

2 

2A(u,v) =T(v,u) = (uAv)Q(x). (140) 

n — 1 

Consider a bundle map K : M x R m — > TM which gives rise to the metric 
compatible connection V (theorem 11). Consider the s.d.e. 

dx T = K{X T ) odW T . (141) 

Its generator is (c.f. [27]) 

1 1 m 

-trace(V 9 ) 2 + - ]T V 9 iT(iT) (142) 

2 2 i=i 



while by assumption we have 

1 1 m 1 

#o(<7, Q) = -trace(V^) 2 + (- ]T V^V* - Z) = -trace(V 9 ) 2 + Q (143) 

i=i 

The required result follows after we show 

53 V 9 iT(iT) = - 53 D(K\ K l ) = traced- -) (144) 
i=l i=l 

equals 2Q. For this we note that for all v £ TM, 

m m m 

g(£D(K i ,K i ),v))=g(£S(K i ,K i ),v)) = - ]T g(f(v, K% IT) = 2g(Q,v) 

i=l i=l i=l 

(145) 

Consequently 

trace(V 9 ) 2 + - ]T V 9 (iT)(lT) = trace(V 9 ) 2 + Q = H (g, Q) (146) 
2 i=i 

and the if so constructed is the required map.c.q.d. 

Remarks 8: Of course in the above construction, it is innecessary to start 
with an arbitrary metric-compatible Cartan connection, only the trace-torsion as 
proved already matters. 
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12.4 Navier-Stokes Equations Is Purely Diffusive In Any Dimen- 
sion Other Than 1 

We recall that for any dimension other than 1, the Navier-Stokes equation is a 
diffusion process which arises from a RCW connection of the form (34) with metric 
given by 2i/g, where g is the original metric defined on TM, and torsion restricted 
to the trace-torsion given by ^u T , where r > 0: we shall call this connection 
the Navier-Stokes connection with parameter v, which we shall denote as \/ NS '^. 
Let us then consider a bundle map K T : M x R m — > TM, for r > for such a 
connection; from theorem 11 we know it exists. Therefore, from Theorem 12 we 
conclude that: 

Theorem 13: For any dimension other than 1, the random lagrangian rep- 
resentations given by (68) admit representations as a Stratonovich s.d.e. without 
drift — u T term 

dx T = K t (x(t)) o dW T . (147) 

Remarks 9: Thus, we have gauged out the velocity in the dynamical repre- 
sentation for the fluid particles. Of course, the new diffusion tensor K T (r > 0) de- 
pends implicitly in the velocity of the fluid as well as in the kinematical viscosity. 
Indeed, K can be in principle computed from the knowledge of the Navier-Stokes 
connection with parameter v, by solving the equation 

v ns-^ z = K(x)d(K\.)Z(.)(v),v € T X M, Z € F(TM). (148) 

Remarks 10: This last theorem deserves further examination. As well 
known, the Navier-Stokes equations is a mechanism of "competition" between 
the linear diffusion term v/S\u T which at the level of the fluid flow described by 
the diffusion tensor term (2v) ?X and the drift non- linear term PL Ut u t , described 
at the level of the fluid flow by the drift vectorfield — u T . The long-time existance 
of the representations of Navier-Stokes depends on the prevailance of the diffusion 
term, so that the non-linearity which feeds the fluid with energy would fade out. 
The representation given by Theorem 13 assures that this is the case for any 
physically interesting dimension, without distinction between 2-dimensional and 
3-dimensional manifolds, since we can always find a non-linear representation 
for the random fluid flow which incorporates both the diffusion tensor and the 
velocity into a new diffusion tensor! (We can think here in the dispute on pre- 
Socratic Greek philosophy between the followers of Parmenides and Heraclitus, 
and thus name (147) as the Parmenidean representation: in this representation 
the fluid flow does not depend explicitly on the fluid velocity.) 

Remarks 11: Just like in theorem 13, we can represent the lagrangian 
random paths for the kinematic dynamo problem as a purely diffusive process, 
by putting v m instead of v. 
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13 Final Observations: 



The method of integration applied in the previous section is the extension to 
differential forms of the method of integration (the so-called martingale problems) 
of elliptic and parabolic partial differential equations for scalar fields [24,31]. In 
distinction with the Reynolds approach in Fluid Mechanics, which has the feature 
of being non-covariant, in the present approach, the invariance by the group of 
space-diffeomorphims has been the key to integrate the equations, in separating 
covariantly the fluctuations and drift terms and thus setting the integration in 
terms of covariant martingale problems. The role of the RCW connection is 
precisely to yield this separation for the diffusion of scalars and differential forms, 
and thus the role of the differential geometrical structure is essential. 

The solution scheme we have presented gives rise to infinite particle random 
trajectories due to the arbitrariness of the initial point of the Lagrangian paths. 
This continuous infinite particle solution is exact and we have actually computed 
its expression. Actually, to integrate NS we choose a finite set of initial points 
and we take for Slo a linear combination of 2-forms (or area elements in the 2- 
dimensional case) supported in balls centered in these points, the so-called many 
vortices solutions; one can choose the original fo(x) so that J7o is supported in 
these balls and these localizations persists in time. Thus the role of the potential 
term in the Buttke magnetization 1-form in the expression (56) is to push the 
vortices to be confined on predetermined finite radii balls; see Chorin [1]. This 
requires that convergence to a solution of NS be proved in addition. 

A new approach to NS as a (random) dynamical system appears. Given a 
stationary measure for the random diffeomorphic flow of NS given by the sta- 
tionary flow of equation (67), one can construct the state space of this flow and 
further, its random Lyapunov spectra. Consequently, assuming ergodicity of this 
measure, one can conclude that the moment instability of the flow is related to a 
cohomological property of M, namely the existance of non-trivial harmonic one- 
forms 4>, which are preserved by the vectorfield u of class C 2 , i.e. L^cf) = di^cp = 0; 
see page 61 in [27]. We also have the random flows {v T A v T : r > 0} and W 2 ~ u ° 
on TM A TM of Theorem 7 and Theorem 10 respectively, which integrate the 
linear equation for the vorticity. Concerning these flows, the stability theory of 
NS (67) requires an invariant measure on a suitable subspace of TM A TM and 
further, the knowledge of the spectrum of the one-parameter family of linear 
operators depending on u, H2(2ug, —^u T ). The latter may play the role of the 
Schroedinger operators in Ruelle's theory of turbulence, which were introduced 
by linearising NS for the velocity as the starting point for the discussion of the 
instability theory; see article in pages 295 — 310 in Ruelle [48]. 

Finally, it has been established numerically that turbulent fluids resemble 
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the random motion of dislocations [4]. In the differential geometric gauge theory 
of crystal dislocations, the torsion tensor is the dislocation tensor [12], and our 
presentation suggests that this analogy might be established rigorously from the 
perspective presented here. We would like to remark that the results presented 
in this article, are part of a more general program of formulation of gravitation, 
quantum mechanics and irreversible thermodynamics, in terms of stochastic dif- 
ferential geometry, developed by the author. 

14 Appendix 

We shall review some basic concepts of the probabilistic, analytical and geomet- 
rical realms. 

Let {T T ■ r > 0} be a family of sub cr-fields of a cr-field J 7 in a probability 
space (ft, J 7 , P). It is called a filtration of sub cr-fields if it satisfies the following 
three properties: i)JF s C T T if s < r; ii)n e> o^v+ e = T T , and hi) each T T contains 
all null sets of T 

A stochastic process x t ,t G T, with T a time-set, say [0, oo), the interval 
[0, T] or the real line, is called {T T ) — adapted if for each r, x T is jF T -measurable. 
Predictable sets are subsets of [0, oo) x R, which are elements of the smallest 
cr-algebra relative to which all real ^v-adapted, right-continuous processes with 
left-hand limits are measurable in (t,u>). A stochastic process x : [0, oo) i— > S, 
where S is a measurable space with cr-algebra B is called predictable if, for any 
Borel subset B £ S, {(t,uj),x(t,lu) € B} is predictable. 

A positive random variable t is called a stopping time (with respect to the 
filtration {T T : r > 0} if for all < r, {t < r} G T T . This concept is used to 
indicate the ocurrence of some random event. 

The conditional expectation of a real-valued random variable X with 
respecto to a sub-cr algebra Q of T is denoted by E(X\Q). 

Let x T be a (J^ T )-adapted real-valued process such that for each r, x T is 
integrable. It is called a martingale if it satisfies: E[x T \T s } = x s a.s. for any 
t > s. Furthermore, x T is called a local martingale if there exists an increasing 
sequence of stopping times {r n } such that P(r n < T) — > as n — > oo, and 
each stopped time x 1 ^ = x rATn is a martingale, where r A r n = mm{r, r ra }. A 
martingale is obviously a local martingale (set r ra = T, for all n. Finally, x T is 
called a semimartingale if it can be decomposed as the sum of a localmartingale 
and a process of bounded variation. 

This decomposition, which in the course of the presentation of the rules of 
stochastic analysis on manifolds for differential forms appears explicitly in the Ito- 
Elworthy Formula, sets the integration of Navier-Stokes equations and the kine- 
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matic dynamo problem as the solutions of the so-called martingale problems 

(after Stroock and Varadhan [30]). 

Let D be a domain in Euclidean space R d and let R l another Euclidean 
space (eventually d = I). Let m G N; denote C m = C m (D;R l ) the set of all 
maps f : D —> R l which are m-times continuously differentiable. For the multi- 
index a = (cci, . . . , a,i) G N d , define 

D °=(5?)^^' with|a| = g^ (I49) 

Let K be a subset of D. Set 

\\f\\m,K = su PxeK J( X \ + Yl su Px&K \D a f(x)\. (150) 

+ l<|Q|<m 

Then C m (D;R l ) is a Frechet space by seminorms {||.|| m; jr : K are compacts in 
D}. When K = D we write ||.|| m; x as Now let 5 such that < 5 < 1. 

Denote by C m;<5 = C m < s {D; R l ) the'set of all / G C m such that L> a /, \a\ = m are 
5-Holder continuous. By the seminorms 

imi imi \- D a f(x)-D a f(y) 

\\f\\m+5;K = \\f\\m;K + 2_, su Px^K,x^y , _ ,g , (151) 

a=rn \ X V\ 

is a Frechet space, the so-called space of 5 — Holder continuous C m mappings. 

When D = K we write ||.||m-i-<5;A~ ^ \\'\\m-\-6' Denote further as ' the set 
{/ G C m ' 5 : ||/|| m +5 < oo}. A continuous mapping f(r,x), x G D,t G T} is 
said to belong to the class C m;<5 if for every r, /(r) = /(r, .) belongs to C m;i5 and 
||/ (r^lm+^i^ is integrable on T with respect to r in any compact subset K. If 
the set K is replaced by D, / is said to belong to the class C™'" 5 . Furthermore, 
if ||/(r)|| m+ 5 is bounded in r, it is said to belong to the class C^" 5 . 

Consider the canonical Wiener space £1 of continuous maps uo : R —>■ 
R d ,iu(0) = 0, with the canonical realization of the Wiener process W(t)(lu) = 
cj(t). Let (f) SjT (x,u>), s,t G T, x G R d be a continous ^-valued random field 
defined on (fi, J 7 , P). Then, for almost all uj G f2, (f) S;T (co) = (f) S;T (.,u) defines a 
continuous map from R d into itself, for any s, r. Then, let us assume that 4> StT 
satisfies the conditions: (i) (j) S)U {io) = (j) T)U {u))(j) StT {uj), holds for all s,t,u, where in 
the r.h.s. of (i) we have the composition of maps, (ii) (f) s , s (Lo) = id, for all s, where 
id denotes the identity map, (hi) <^> SjT (w) : R d — > R d , is an onto homeomorphism 
for all s,t, and (iv) (f) S;T (x,uj), is an onto homeomorphism with respect to x for 
all s,t, and the derivatives are continuous in (s,t,x). 
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Let SiT (a;) _1 be the inverse map of (f> SjT (u). Then i) and ii) imply that 
<p T<s (uj) = SjT (t<j)~ 1 . This fact and condition hi) show that SiT (u;) -1 is also 
continuous in (s,t,x) and further condition iv) implies that (p SjT (u;)^ 1 (x) is k- 
times differentiable with respect to x. Hence (p s ^(uj) : R d — > R d is actually a 
C fc -diffeomorphism of M, for all s, r, the so-called random diffeomorphic flow. 
We can regard 4> StT (iv)~ 1 (x) as a random field with parameter (s,t,x) which is 
often denoted as </>~*(x,u;). Therefore, 

= <j> T ,s(x) (152) 

holds for all s, r, x a.s.When we choose the initial time s of (j> SjT (x, to be 0, we 
shall write <J) t (x,lo). 
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